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18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting # 8: Green Functions 



1. Green functions for domains Q 

Our goal in this section is to derive an integral representation formula for the solution to Poisson's 
equation on domains fl c M". Specifically, we will study the boundary value Poisson PDE 



(1.0.1) 



Au{x) = f{x), 
u(x) = g{x), 
We first state a basic existence theorem. 
Theorem 1.1 (Basic existence theorem). Let g be a bounded Lipschitz domain, and let g e 



C{dil). Then the PDE (1.0.1) has a unique solution u e C^(fi) nC{^l). 



Proof. This proof is a bit beyond this course. □ 

We now define the basic object that will play the role of a fundamental solution on a domain Q. 

Definition 1.0.1. Let Q c M" be a domain. A Green function in Q is defined to be a function of 
{x,y) eQ xQ verifying the following conditions for each fixed x : 



(1.0.2) AyG{x,y) = 6M''=Hy-^), 

(1.0.3) G{x,a) = 0, 

Let's now connect G(x,y) to $(x-?/). 

Proposition 1.0.1. Let $ be the fundamental solution for A in 
Then the Green function G{x,y) for Q can be decomposed as 

(1.0.4) G{x,y) = <^{x-y)-(P{x,y), 

where for each x eQ, (f)(x,y) solves the Dirichlet problem 



y e ^, 

and let £ be a domain. 



(1.0.5) 
(1.0.6) 



Ay(j){x,y) = 0, 

a) = $(x - cr), 



y e ^, 



Proof. As we have previously discussed, A$ = 6. Also using (1.0.5), we compute that 



(1.0.7) Ay{<!>{x-y)-(P{x,y))=Ay^x-y)-Ay(l){x,y) = 5{y-x). 

Therefore, $(x - y) - (p{x,y) verifies equation ( |1.0.2 ). 
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Furthermore, using (1.0.6), we have that <l>(x - a) - (j){x,a) = whenever a e dfl. Thus, $(a; 
y) - 4>{x,y) also verifies the boundary condition (1.0.3). 



□ 

The following technical proposition will play later in this section when we derive representation 



formulas for solutions to (1.0.1) in terms of Green functions. 



Proposition 1.0.2 (Representation formula for u). Let $ be the fundamental solution for A 
in M", and let Q c M" be a domain. Assume that u e C'^{Q). Then for every x e Q, we have the 
following representation formula for u{x) ■ 

(1.0.8) u{x)= f ^{x-y)Ayu{y)d''y- f ^{x - a)V fj^^Ma) da + f u{a)V f,^^)^{x - a) da . 

" , ' , ' 

single layer potential double layer potential 

Proof. We'll do the proof for n = 3, in which case = -4^- We will also make use of Green's 

identity. Let B^(x) be a ball of radius e centered at x, and let fi^ '^=^ ^l\B^(x). Note that dil^ = 
dfl u -dB^{x). Using Green's identity, we compute that 

(1.0.9) / -/^u{y)d'y= / -V fju{a) - u{a)V J-. -)da 

Jn,_\x-y\ Jdii^\x-a\ ^\x - a\' 

= / 1 j^u(a') da - \ u(a)"^ A-. -Ada 

JdVL\x-a\ Jan \\x-a\^ 

-\/j;>u{a)da+ I u{a)\I A- -Ada. 

dB^{x)\x-a\ JdB,(x) ^\x-a\' 

In the last two integrals above, N{a) denotes the radially outward unit normal to the boundary of 
the ball B^{x). This corresponds to the "opposite" choice of normal that appears in the standard 
formulation of Green's identity, but we have compensated by adjusting the signs in front of the 
integrals. 



Let's symbolically write (1.0.9) as 
(1.0.10) L = R1 + R2 + R3 + m. 

Our goal is to show that as e i 0, the following limits are achieved: 

• L^-4n $(a; - y)Ayu{y) d^y 

• Rl Ait X single layer potential 

• R2 -Att X double layer potential 

• i?3 ^0 

• i?4 -47rn(x). 

Once we have calculated the above limits, ( |1.0.8 ) then follows from simple algebraic rearranging. 
We first address L. Let M = maXy^-^ Au{y) . We then estimate 



:i.o.ir 



[ ^^Auiy)d^y- [ -^Auiy)d^y < [ -^\Au{y)\d^y 
jQ\x-y\ Jn^\x-y\ JB,{x)\x-y\ 



<M 



'b,(x) \x - y\ 



d^y ^ as e i 0. 



MATH 18.152 COURSE NOTES - CLASS MEETING # 8 



3 



This shows that L converges to |^3^Am(|/) (t^'y as e j 0. 

The hmits for Rl and R2 are obvious since these terms do not depend on e. 

We now address i?3. To this end, Let M' = max^g^^ |Vii(y)|. We then estimate R3 by 



(1.0.12) 



\x - a 



da 



< f -M'da= 47re^ xe-^M' ^ as e 4 0. 

JdB^x) e ' — . — ' 

surface area of dBe(x) 



We now address R4. Using spherical coordinates (r,9,(f)) e [0, oo) x [O.tt) x [0,27r) centered at 
X, wc have that da = sin 9 dO d(f). Therefore, fg^^f^^-^ \x-a\'^ ^"^ ~ f(pe[o 27r] fee[o it] ^dOdcf) - Att. We now 
estimate 



;i.0.13) 



—R4:-\-u(x)] = u(x) + — / u(a)V Mt^A ] ndcr 

/ {u{x) - u{a))[ -) da 

JdB^ix) ^ '^\x-aY'f 

< — / \u{x) - u{a)\i- rAda 

Att JdBeix) v|a;-(7p/ 

X X 
<— max \u(x)-u(a)\ / \ ;7;]da 

47r cre9B.(x) ' JdB,(x)\\x - a\^ ^ 



1 

An 



< max \u(x) - u(a)\ -y as e I 0. 



This shows that RA -Anu{x) as e J, 0. 



□ 
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